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In  cctaputing  the  rate  of  advance  of  the  crystallization  front  during 
solidification  of  a nsetal  ingot,  one  usually  assumes  that  tlie  themal  con- 
ductivity X ^d  the  heat  capacity  c of  a solid  metal  are  independent  of  the 
temperature.  Die  average  valuea  of  the  thernophysical  parmnetera  of  the 
substances  X and  ~c  in  the  given  temperature  range  are  used  in  the  computa- 
tion. As  the  range  is  extended,  the  actual  values  of  X end  c may,  at  certain 
temperatures  T,  differ  considerably  from  X and  c.  The  error  thus  introduced 
into  the  result  retnaina  undetermined,  and  therefore  the  validity  of  using 
X and  X instead  of  the  actual  thermophyaical  paramctera  of  the  metal  is 
doubtful.  To  solve  the  problem  of  the  possibility  of  coBsputing  the  solidifi- 
cation r?te  of  the  metal  with  sufficient  accuracy  while  limiting  ouraelves 
to  ^ and  c,  we  must  have  a more  accurate  method  enabling  us  to  allow  for 
the  temperature  dependence  of  X and  c. 

In  this  paper  we  try  to  develop  such  a method,  ahich  could  be  used  for 
various  dependences  of  X and  c on  T.  For  concreteness  we  assume 

X = Xj-fp7',  c = c^  + o)T,  (1) 

The  relation  between  X,  c,  and  T corresponding  to  equations  (1)  applies 
to  many  metals  and  alloys  in  a wide  temperature  range.  For  most  metals,  we 
ha»»  p<0,  o)>0. 


We  shall  assume  that  there  is  no  superheating  cf  ths  sslt  and  that  the 
temperature  of  the  metal  depends  only  on  the  distance  from  the  surface  of 
the  ingot  and  on  the  tiina  t (reference  2).  The  equation  describing  ''he 
variation  of  the  temperature  in  the  solidified  crust  of  the  metal  has  the 
form 


^ + pT)  g]  = 7 (c, o>  D : 


iriiere  y density  of  the  metal; 


T(X.  0)  = Tcr  tor  x>0,  (3) 

where  is  the  crystalliiation  temperature  of  the  metal. 

If  a constant  temperatnre  is  preserved  at  the  surface  of  the  ingot  dur- 
ing crystallization,  we  have,  for  t > 0, 

T{0.  0 = To  , (4) 
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where  To  is  the  temperature  of  the  surface  of  the' ingot. 

Oi  the  crystallization  front  * » y(t),  where  y(t)  is  the  coordinate  of 
the  crystallization  front,  increasing  with  time, 

T[y{t),  (5) 

From  the  condition  of  heat  balance  at  the  melt-liqpiid  [sic]  phase  boundary, 
we  obtain 


dy 

da 


/dT\ 

Qy 


(-) 

\dx  Jj,=.y{t)  ’ 


(6) 


«dtere  Q is  the  heat  of  crystallization  per  gram. 
We  introduce  the  quantities: 


f = Cu  + u»  ■ 


a = 


0 = 

2 

/ J ^0  + ^cr  \ . 

1 

u 

X 

— I 

TC 

0t|  — _ 

2X 

io). 

> 

<>>  (Ter  — T^) 

fits  * j 

2c 

20 


Ihe  conditions  of  our  problem  can  be  written  as 


d 

dx 


We  shall  seek 


[(1  (1 0(0,  0 = --i; 

at  g \ nX 


(7) 


(8) 


0 {X,  t) 0 (a);  a = = 2at, 


where  3 is  an  unknown  constant. 
Hence 


(9) 


de  1 , 

(l+«ie)^  = 0; 

(10) 

0(0)  = 

— 1; 

(11) 

A (H)  = 

= 1; 

_ ^ «l/ 

rfe\ 

(12) 

? 1 

du  /u=9  ■ 

(13) 

The  case  ■ 0 corresponds  to  taking  the  average  values  for  X 

and  c in  the  entire  temperature  range  from  T to  Tq.  Then 


^ out 

(14) 

e(0)  = -l; 

(15) 

0(P)  = 1; 

(16) 

(17) 
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here  0 is  the  value  of  0 conputed  from  the  average  values  of  X and  c. 

The  solution  of  Eq.  (14)  with  due  regard  for  the  conditions  in  (15)  and 
(16)  has  the  fom 

2ert  (u  / V2)  . 2 r 

erf(P/1^2)  (18) 


0^ 


Condition  (17)  results  in  the  transcendental  equation  determining  the 
value  of  0; 


2VT 


=— ?-g-(3)*/*  y 

' iFrj 

W2J 

IT. 

1.3  o-7.-.(2*41)  ' 

Since 


it  follows  from  (19)  that 

f = + r43(P)‘ + 1^:7 (?)*  + ••  • 

The  solution  of  Eq.  (10)  can  be  written  in  the  form  of  the  series 


e(u)==  M^nj  — 


It  can  be  seen  from  conditions  (12)  and  (13)  that 


0( 


=--1;  ( 


f[e  \ ^ «?P 

au  1 -h  «i ' 


(19) 


(20) 


(21) 


(22) 


(23) 


By  successive  differentiation  of  Eq.  (10),  we  find  the  derivatives  of 
the  function  6 (u)  of  higher  order  for  u “ 0: 

r_L±i^L_.  a 4. 

\ du«  ^ [ (1  + «i)‘  " (1  + ®i)*  J ’ 

4-  B»r  4a.  + 3a.«,-(x,  , 3«» 

_ M f3(t+«dV.  4.  K+4..g,-2«,  _ jn  -H^)« 

U«V»-9  “ P L (1 ‘ (l+«d*  + 


(1  + + 7«|04  — 4a,) 

(1  + «i)‘ 


+ 

etc. 

According  to  condition  (il) 


4-  «i(25«i  + 18a.«,  — 7a,)  . , -i 

V (!+«,)•  ^ ^ (1  + «i)’  ^ J 

1) 

«=o  ' ““  '«— 8 


(26) 


(27) 


Substituting  equations  (24),  (25),  (26),  etc.,  into  (27),  we  obtain 
after  certain  transfomations  a series  of  even  powers  of  the  reqpiired 
constant  0: 

Y (1  + «i)  = ?*  + ?*  [ 3 (1  4 «,)  + 2 (1  4 «.)•  + P*  [15  (1  4%* 

• 5«i  + 4xi04  — «,  , ,1  , cf  (l+«i)*  I (I4a«)(31»i422aia,— , 

12(1  4 a.)'  ^^2(14«.)‘^  J‘^*^1105<14-.)''^  180  (1  4 ai)‘  ^ 


g,(19«,4  ISaig,  — 4a,)  , ^a*  -1 

^ 30  (1  4 ^ 8 (1  + a,)*  V j + 


(28) 
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It  is  essy  to  see  thst  for  OC],  * * 0 the  scries  (28)  is  transfonied 

into  the  series  (21). 

When  the  vslnes  <Xi  and  are  giren  concretely,  the  conrergence  of  the 
series  (28)  is  studied  by  the  usual  methods. 

Expression  (28)  is  a transcendental  equation  from  whxch  it  is  possible 
to  determine  the  rslue  of  6 graf^ically  or  by  the  method  of  inrersion  of 
the  poirer  series.  When  3 is  knosm  it  is  easy  to  calculate  the  rate  of  ad- 
rance  of  the  crystallization  front  during  the  solidification  of  the  metal 
ingot: 

dy  ^ 

V2i  * (29) 


In  the  case  of  a more  complex  dependence  of  X and  c on  T,  the  scheme  of 
the  solution  of  the  problem  remains  unchanged,  except  that  the  cocnputation 
formulas  will  be  different. 

There  may  be  cases  when  the  computation  can  be  made  according  to  (21) 
with  great  accuracy  and  the  value  of  3 is  practically  the  same  as  3* 

Let  us  calculate  the  crystallization  rate  of  an  aluminum  ingot.  In  this 
caae^  ws  have  Q “ 95.5  cal/gm,  “ 659°  C,  Y * 2.7  gm/cm**.  The  de- 
pendence of  X and  T for  alumintim  hsa  been  reported.  Analytically, 

the  dependence  of  \ and  c on  T in  a wide  temperature  range  can  be  de- 
scribed by  the  formulas 


cal 


c=(0, 2096  + 0, 0001 TC) 


cal 


tP>  X deg 


X~(0.5-0.00025TX)^^.^^j.^, 

If  we  aastiae  that  Tq  * 100°C,  we  have  Oti  = -0,1725;  Otj  " 0,1129;  and 
q ■ 1.3803.  From  the  series  (28K  »e  obtain  6 * 1.011.  For  ttj  * flj  ■ 0, 
the  computation  gives  ^ * 0.989. 

In  the  case  under  study  3/0  ■■  1.023;  i.e.,  the  error  in  determining  3 
caused  by  the  temperature  dependence  of  the  theimophysical  parameters  of 
aluminum  is  somewhat  larger  than  2 percent  of  the  computed  quantity.  Thus, 
in  calculating  the  solidification  of  the  aluminum  ingot,  we  can  assusm 
with  sufficient  accuracy  that  3 ^ 

There  are  possible  cases  in  the  crystallization  of  metals  (particularly 
metals  having  s hi^  melting  point)  when  the  difference  between  3 and  0 
is  much  more  substantial  than  in  the  above  example. 
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